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Abstract. We consider gravity duals to d+1 dimensional quantum critical points with 
anisotropic scaling. The primary motivation comes from strongly correlated electron systems in 
condensed matter theory but the main focus of the present paper is on the gravity models 
in their own right. Physics at finite temperature and fixed charge density is described in 
terms of charged black branes. Some exact solutions are known and can be used to obtain 
a maximally extended spacetime geometry, which has a null curvature singularity inside a 
single non-degenerate horizon, but generic black brane solutions in the model can only be 
obtained numerically. Charged matter gives rise to black branes with hair that are dual to the 
superconducting phase of a holographic superconductor. Our numerical results indicate that 
holographic superconductors with anisotropic scaling have vanishing zero temperature entropy 
when the back reaction of the hair on the brane geometry is taken into account. 



1. Introduction 

We consider d + 2 dimensional gravity models that are dual to quantum critical points with 
anisotropic scaling in d + 1 dimensions, 

t -> X z t, x ->• Ax, (1) 

with z > 1 and x = (xi, . . . , Xd)- The d + 2 dimensional bulk action consists of three parts, 

& = "^Einstein-Maxwell + "^Lifshitz + "^matter ■ ( 2 ) 

The first term is the standard action of Einstein-Maxwell gravity with a negative cosmological 
constant, 

•^Einstein-Maxwell = / d"+ 2 xj=g (r~2A- ~ V*") • (3) 
This is followed by a term involving a massive vector field, 

%fshitz = - / d d+2 x^-g (j V" + , (4) 

whose sole purpose is to provide backgrounds with anisotropic scaling. This is a modified version 
of the holographic model of pQ , which was formulated in four-dimensional spacetime and obtained 



anisotropic scaling by including a pair of coupled two- and three- form field strengths Q Upon 
integrating out the three-form field strength, the remaining two-form becomes a field strength of 
a massive vector [3] and in this form the model is easily extended to general dimensions. Finite 
temperature corresponds to having a black hole in the higher-dimensional spacetime and with 
the Maxwell gauge field, added in [3], the black hole can carry electric charge, which is dual to 
a finite charge density in the lower-dimensional theory. Finally, we consider matter in the form 
of a scalar field, which is charged under the Maxwell gauge field but does not couple directly to 
the auxiliary massive vector field, 

Matter = ~\ J d d+2 x^~g ((d^* + iqA^*) (d^ - iqA^) + m 2 0>) . (5) 

At low temperature there is an instability for charged black holes in the model to grow scalar 
hair, which corresponds to the superconducting phase of holographic superconductors with z > 1 
asymmetric scaling [JJ, as described in Section [3] below. 

The equations of motion obtained from the action ^) consist of the scalar field equation 

(V»-i q A> 1 )(V fi -iqA fi )(f>-m 2 (f> = 0, (6) 

along with the Einstein equations, the Maxwell equations, and field equations for the auxiliary 
massive vector, 

/~i i a _ ^Maxwell , ^Lifshitz , T-imatter tn\ 

V u F vp - = j matter , (8) 
V v jw = c 2 A». (9) 



The asymmetric scaling symmetry (TJJ), sometimes referred to as Lifshitz scaling, is realized in a 
d+2 dimensional spacetime, 

ds 2 = L 2 (-r 2z dt 2 + r 2 dx 2 + ^ , (10) 
whose metric is invariant under the transformation 

r 

t -> X z t, x ->• Xx, r -t — . (11) 

A 

Length dimensions are carried by the characteristic length L while the coordinates (t, r, x) are 
dimensionless. The scale invariant Lifshitz geometry (|10p is a solution of the equations of motion 
when L is related to the cosmological constant A via 

z 2 + (d-l)z + d 2 



and the mass of the auxiliary vector field is fine-tuned to c = yzd/L. In this solution the 
Maxwell field vanishes, = 0, but the massive vector field has the background value 



At = \j^-^L r z , A Xi =A r = 0. (13) 



See also [2] for early work on gravitational backgrounds with anisotropic scaling. 



2. Charged black branes 

In order to study finite temperature effects in the dual strongly coupled field theory we look 
for static black brane solutions of the equations of motion (JSJ) - ([9j) which are asymptotic to the 
Lifshitz fixed point solution given by (|10p and (|13p . From now on we set L = 1 and consider a 
metric of the form 

ds 2 = _ r 2z f{r) 2 dt 2 + r 2 dS 2 + 9(T? dr 2 _ (M) 

An asymptotically Lifshitz black brane with a non-degenerate horizon has a simple zero of both 
/(r) 2 and g(r)~ 2 at the horizon, which we take to be at r = ro, and f(r), g(r) — >• 1 as r — >■ oo. 
It is straightforward to generalize this ansatz to include black holes with a spherical horizon or 
topological black holes with a hyperbolic horizon but it is the flat horizon case (|14j) that is of 
direct interest for the gravitational dual description of strongly coupled d+1 dimensional field 
theories. In a static electrically charged black brane background the Maxwell gauge field and 
the massive vector can be taken to be of the form 



A, = r z f(r){a(r),0, . . . , 0} ; A» = ]J ^^r z f(r){a(r), 0, . . . , 0} . (15) 
with a(r) — > and a(r) — >• 1 as r — > oo. 



2.1. Field equations for static configurations 

For static configurations the equations of motion can be expressed as a first order system of 
ordinary differential equations. Introducing a scale invariant radial variable u = log(r/ro) and 
writing ' = jfjr, the scalar field equation becomes 



(d+z) \ + (m 2 



2 2\ 2 

J a ) g c 



while the Maxwell equations and the equations of motion for the massive vector reduce to, 



a + —a 
a + J a 



-za + g/3, 
-dp + q 2 g <j) 2 a , 
-z a + z g b , 
-db + dg a . 



(16) 
(17) 

J 

(18) 
(19) 
(20) 
(21) 



The functions Xi fti an d b are defined via (fT6|) . (fT8)) . and ([20]) . respectively. The Einstein 
equations can also be written in first order form, 



9 / 
9 f 
f 

2dL 



9 2 a 2 - 1 ) + 



( 2 , 2 2 2-2 



(22) 



d (\-d-2z) + ^ + g 2 



2 



z—1) (^da 

+- (q 2 a 2 - m 2 ) (f> 2 + z 2 + + d 



(23) 



The field equations (|16p - ()23p are manifestly invariant under the scaling (jlip and the Lifshitz 
fixed point solution is given by/ = g = a = 6 = l and a = f3 = (p = x = 0. 



In the absence of charged matter, the Maxwell equation (|19p integrates to j3(u) = pe~ du . 
The integration constant p = p/r$ is proportional to the electric charge per unit d- volume of the 
black brane, which in turn corresponds to the charge density in the dual field theoretic system. 
For given d and z > 1, the remaining field equations then have a one parameter family of black 
brane solutions, labelled by p. A neutral black brane without scalar hair has p = while the 
extremal limit is given by p -»■ ±y/2(z 2 + (d-\)z + <£ 2 ). A non- vanishing charged scalar field 
changes this picture as discussed below. 

Equation (JTHJ) can easily be solved in the Lifshitz geometry (fTUj) without matter, giving 



1 7. „—du 



pe- au if z^d, 



a(u) = l P ° . ^1^)7 " (24) 

v ; \ ~pe~ du + pue~ du ii z = d, y J 

where the integration constant p = p/r^ corresponds to having non- vanishing chemical potential 
in the dual system. In general, we are not working with the Lifshitz background but with 
solutions that are only asymptotically Lifshitz as u — > oo. However, as long as the value of z 
isn't too highjl the asymptotic behavior of the gauge potential carries over from the Lifshitz 
background to the more general case, and one can read off the charge density and chemical 
potential in the dual field theory from the leading two terms in the expansion of a at large it. 
Calculations in this paper refer to fixed p, corresponding to a fixed density of charge carriers in 
the dual system, but one can also work at fixed chemical potential. 



2.2. Numerical solutions 

Black brane solutions at generic z > 1 can be obtained using numerical techniques similar to 
those of [jJJ. The field equations are integrated numerically starting near the black hole, with 
suitable initial conditions and proceeding out towards the asymptotic region. For a regular, 
non-degenerate horizon we require the functions that appear in the metric ansatz (|14p to behave 
as 

f{u) = Vu(J + fiu + ...), g(u) = — ( go + 9lU + ...), (25) 

y/U 

near u = 0. The appropriate near-horizon behavior of the remaining field variables can be 
worked out and then inserted into the equations of motion to generate initial value data for the 
numerical integration. For any given values of d and z, we obtain a three parameter family of 
initial values, where for instance (fi(0), (3(0), and 6(0) can be taken as the independent parameters. 
The condition that the metric and massive vector approach the Lifshitz fixed point solution (jlOp 
and (fTBj) sufficiently rapidly as u — > oo restricts the solutions further O E]. As a result, 
6(0) is fixed for given 0(0) and /3(0) and one has a two-parameter family of solutions. This 
means in particular that, in the absence of scalar hair, there is a unique (up to overall scale) 
asymptotically Lifshitz charged black brane solution for given d, z and p. 

The Hawking temperature is determined by the near-horizon behavior of the black brane 
metric, 

T H = ^^. (26) 

The full numerical solution of the field equations is required, however, to relate the coefficients 
/o and go to the charge density p and other physical variables of the dual field theory. 



2 At z > 3d non-linear effects give rise to additional terms in ()24|l with a falloff in between that of the charge 
density and chemical potential terms. 



2.3. Conserved charge under radial evolution 

A conserved charge under radial evolution was found in [6] for electrically neutral black branes 
with d = 2 and arbitrary dynamical critical exponent z. Such a conserved charge is useful for 
matching solutions across the bulk geometry from the near-horizon region to the asymptotic large 
u region and also provides a check on numerical solutions. The charge found in [6] generalizes 
to charged black branes with scalar hair in general spatial dimensions, 



D = r ^+d e ^+d) Uf 



+9 



— (X 2 - 2d(d+l)) - 2d(z-l)ab -da/3 

1 



(z-l)(da 2 -zb 2 ) + z 2 + {d-l)z+d 2 - -(/3 2 + mV - q 2 a 2 (p 2 ) 



(27) 



It is straightforward to check that J^-Do = when the field equations (|16p - (j23[) are satisfied. 
The conserved charge is related to thermodynamic state variables of the dual system in a simple 
way. Inserting a perturbative near-horizon expansion of the fields, one finds 



2r, 



d+z 



fo 

go 



32ttST, 



(28) 



where T is the temperature (|26p and S = Tq /4 is the Bekenstein-Hawking entropy density, which 
are identified with the temperature and entropy of the dual field theory. 



2. 4- Exact solution 

As always, it is useful to have explicit analytic solutions to work with. Although Lifshitz black 
branes at z > 1 can in general only be obtained numerically, it turns out that for each value of 
d an isolated z = 2d exact solution can be found gUBji 



1 



a 2 = l-e- 2du , fa = ±V2e- du (l-e~ du ), 



(29) 



with (f) = and p = ±y/2d. It is straightforward to continue the exact black brane metric inside 
the horizon and obtain the globally extended geometry [3]. Define a tortoise coordinate by 



2dr\ d 



log 1 



(l - e~ 2du ) 



(30) 



and then transform the (t, u) variables to a pair of null coordinates 



U 



V = exp drl d {u*+t) 

In the new coordinate system the metric is given by 

-dUdV 



cxp 



drl d (u*-t) 



ds z 



+ 



Tq dx 2 



d 2 {i + uvy (i + uvy/ d ' 



(31) 



(32) 



and is manifestly non-singular at the horizon, which is located at UV = 0. There is a null 
curvature singularity at UV — > oo, which corresponds to r — > in the original coordinate 
system. The asymptotic region r — > oo corresponds to UV —)•—!. By a further transformation 



3 The corresponding exact solutions for d — 2, z — 4 black holes with a spherical horizon and topological black 
holes with a hyperbolic horizon were also found in [4]. In the limit of vanishing electric charge these black hole 
solutions reduce to the previously discovered z — 4 black hole solution of [6]. 



Figure 1. Global geometry of an asymptotically 
Lifshitz black hole. The null singularity at r = 
is indicated by the jagged lines. 



to new null variables P, Q, defined through V = tan 2 ^, U = tan 2^, the global geometry can 
be represented by a simple diagram shown in Figure [TJ Each point in the diagram represents 
an entire G?-volume parametrized by x. The global diagram in Figure [T] differs from standard 
Carter-Penrose conformal diagrams in that the boundary at r — > oo is not conformally flat. 
This is a consequence of the scaling asymmetry between t and x and is readily apparent in the 
globally extended metric (f3"2j) . 

When z = 1 the auxiliary massive vector field An can be consistently set to zero and the 
field equations then reduce to those of Einstein-Maxwell gravity with a negative cosmological 
constant. In this case, there is a well known exact solution, the AdS-Reissner-Nordstrm black 
brane, which has a timelike curvature singularity inside an inner and an outer horizon. The 
interior geometry of the exact z = 2d black brane is markedly different with a null curvature 
singularity at r = and no smooth inner horizon. 



3. Holographic superconductors with asymmetric scaling 

We now consider charged black branes with hair. Static spherically symmetric solutions of the 
scalar field equation ([6]) have the asymptotic form <f>(u) — > c_(e _A_n + ...) + c + (e~ A+u + ...), 
with 

A± .^£ ± ^£) 2 + m2 , (33) 

while the asymptotic behavior of the electromagnetic field field strength is 

(3(u) « ^e~ du + .... (34) 

Working at fixed charge density in the dual field theory, we read the radial location of the horizon 
off from the asymptotic behavior of f3(u) and then the temperature can be obtained from the 
numerical solution for f{u) and g(u) using (|26p . 

In the following we set the scalar mass squared to m 2 = | — {^-^j , which is inside the 
range where there is a choice of two boundary theories [9]. This choice leads to convenient 
values, A± = ± ^, for the dimensions of the operators 0± that are dual to the scalar field 
in each of the two boundary theories. Non-linear descendants of the leading scalar field modes 
are suppressed by 0(e~ 2u ) at u — > oo and this choice of mass squared ensures that the first 
descendant of falls off faster than ip + . 

In order to study holographic superconductivity, we first select some value for d, z and the 
electric charge q carried by the scalar field and then generate numerical black brane solutions for 



a range of initial values (3(0) and (f)(0). We then investigate the asymptotic large u behavior of 
the scalar hair in the numerical solutions. A superconducting condensate corresponds to either 

c+ = 0, (O_) = c_/0, or c_=0, (0 + )=c+^O, (35) 

depending on which of the two boundary theories is being considered [10| 111] . We look for a 
curve in the (3(0) vs. ip(0) plane of initial values at the horizon, for which the corresponding 
black brane solution has vanishing c + (c_), and tabulate the value of c_ (c+) along this curve. 
The temperature is found from the same numerical solutions via (|26p and (|34p . Figure [2] shows 
a plot of c_ vs. T obtained by this procedure for d = 2, z = 2, and q = 1. The results are 
expressed in terms of dimensionless ratios that are insensitive to the overall scale (set by the 
charge density p, which is held fixed at some finite value throughout). 

These results demonstrate that a superconducting condensate can form in systems with 
anisotropic scaling but we have not touched on a number of interesting topics including the 
electric conductivity and magnetic properties of these holographic superconductors. 




Figure 2. Scalar field condensate in a Figure 3. Entropy as a function of 
holographic superconductor at z = d = 2. temperature in the superconducting phase of 

the holographic superconductor in Figure [2] 



4. Zero temperature entropy vs. scalar hair 

In the absence of charged matter, the charged black branes in our model have an extremal limit 
given by p — > ±y^2(z 2 + (d—l)z + d?). It then follows from p = p/r$ that the radial location of 
the horizon ro has a finite value in the extremal limit for fixed charge density p. This in turn 
means that the entropy density S = Tq /4 remains finite in the zero temperature limit indicating 
a macroscopic groundstate degeneracy. 

This conclusion is radically altered when the system is coupled to a charged scalar field. 
In this case the zero temperature limit is approached in the condensed phase and the black 
holes, that are dual to extreme low temperature states, have scalar hair. It is straightforward 
to keep track of the black brane entropy as the temperature is lowered below the critical value 
for forming the superconducting condensate. The result for the same d = 2, z = 2 holographic 
superconductor as was considered in the previous section is shown in Figure [31 Both the entropy 
density and the temperature are normalized to their values at the onset of condensation, S c and 
T c respectively. 

Although the numerical calculations break down before absolute zero is reached, the numerical 
data strongly suggest that S vanishes in the T — > limit for a generic dynamical critical 
exponent z, which is consistent with a non-degenerate ground state in the dual field theory. The 
corresponding result for conformal systems with z = 1 was established in [12) . 



5. Summary 

We have presented an overview of the construction of charged black brane solutions in 
gravity models that realize the anisotropic scaling symmetry that is characteristic of many 
interesting quantum critical points. The motivation for the study of these models comes 
from condensed matter theory, in particular from two- and three-dimensional systems involving 
strongly correlated electrons. The relevance of gravitational models to real world condensed 
matter systems remains highly speculative, but the gravitational approach continues to produce 
effects that are intriguingly similar to what is seen in experiments. A recent example from our 
own work [T3] involves non-Fermi-liquid behavior in the specific heat of anisotropic black branes 
of the type considered in the present paper, which turns out to be qualitatively similar to the 
measured specific in certain heavy fermion metals near a quantum phase transition |14[ I15j . 
While much of the work on gravitational modeling of strongly coupled field theories to date 
involves asymptotically AdS spacetime and an underlying conformal symmetry, it was crucial 
to the success of this particular application to have a non-trivial dynamical critical exponent 
z > 1. This provides impetus for further study of gravity models with anisotropic scaling. 
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